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Abstract. We characterize the jammed state structure of the random sequential 
adsorption of segments of two different sizes on a line. To this end, we define the 
size ratio as a dimensionless quantity measuring the length of the large segments in 
terms of the smaller ones. We introduce a truncated exponential as an ansatz for 
the probability distribution function of the interparticle distance at the jammed 
state and use it to reckon the first four cumulants of the probability distribution 
function. The sole free parameter present in the various analytical expressions 
is tied to the mean interparticle distance from Monte Carlo simulations, while 
the remaining three cumulants are computed without any free parameter and 
compared to Monte Carlo results. We find that the proposed ansatz provides 
results in good qualitative agreement with Monte Carlo ones. 
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Introduction 

The random sequential adsorption model has been used as a paradigm to study 
irreversible phenomena like deposition and adsorption of particles onto a substrate in 
systems where adsorbed particles do not significantly diffuse within the experimental 
time scale [H [51 |3l HJ [5l [6j [7] . Nonetheless, the random sequential adsorption model 
has also been applied to the modeling of apparently dissimilar systems like the cases 
of the classical works on intramolecular reactions taking place in polymers [5J |H] and 
the famous car parking problem [101 111) . The random sequential adsorption model 
has been also applied to study the irreversible adsorption of colloids [T21[T3] and, more 
generally, the deposition of extended objects [UJ [TSl Uni [IZl H] • Recently, interest has 
shifted towards the study of the effect of adsorbing particles of different sizes, with 
research becoming focused on the study of binary mixtures of particles undergoing 
either competitive adsorption or preadsorption. Regarding the latter case, particles 
adsorb during two subsequent stages with each of the particle species attempting 
adsorption at a time (THl [201 [211 [221 [22]) while during competitive adsorption both 
particle species attempt adsorption simultaneously [T7l[T8[|2Ul[24[|25l[26[[27l[28] . 

In general, several physical properties are related to the interparticle distance, 
so proper characterization of the corresponding probability distribution function is 
fundamental. The interest stems from the theoretical perspective in gaining insight 
on the functional dependence of the interparticle distribution even when an analytical 
solution is not available. From the experimental perspective, the availability of the 
interparticle distribution function provides for controlled analysis. In the present work 
we focus on the competitive adsorption of segments of two different sizes on a line and 
characterize various interparticle distribution functions with moments up to the fourth 
order. We propose a novel approach to characterize the gap-size distribution function^ 
based on the systematic use of their cumulants. Comparing the semi-analytical results 
using an ansatz with recently reported Monte Carlo simulation results [28], we find 
good qualitative agreement. 

The paper is organized as follows: In the next section, we start introducing the 
motivation, some terminology, and theory. In Sec. [21 wc introduce the ansatz and 
we compute the various cumulants. We present our results and final remarks in 
Sees. [31 and [H 

1. Motivation and theory 

We consider the competitive adsorption of two size segments on the line, namely 
short segments, which we denote as A segments, and long ones as B. The model 
corresponds to an extension of the continuum random sequential adsorption model 
[21 [3l m [25l 1211 [29], where a segment successfully adsorbs if it does not overlap 
a previously adsorbed one, therefore mimicking an excluded volume, short-range 
interaction. The model has been studied both analytically [251 [27] and by Monte 
Carlo simulations [51 [23 [2H]- Without loss of generality, we rescale the size of the 
system by the length of short segments, so that the length of the smaller segments 
becomes unity while that of the longer segments is R. We term as the size ratio the 
quotient of the lengths of the larger to the smaller segments, which equals R (the length 
of the larger segments in dimensionless units). In general, one can have both segment 



I This is a term that will be interchangeably used to represent the interparticle distribution function. 
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types to arrive at different rates per unit length and per unit time, which represents 
the incident particle flux in the one-dimensional case. In the present work, we take 
the particular case of equal incident fluxes for both types of segments, i.e., on average, 
equal number of segments attempt deposition per unit time and per unit length. Since 
adsorption of segments is irreversible, i.e., in physical systems where surface mobility 
and detachment upon adsorption can be safely neglected, it leads to a jammed state, 
where no segment can fit into any of the available interstitial spaces. The present model 
has been studied in the literature in various limiting cases. For example, for i? = 1 it 
corresponds to the Renyi or car problem result [Tl[llllll[51ll[71[Hl[ni[Tni[ni[lSl[171[30]- 
The case of i? — >■ oo has been previously discussed [28]. The limit of i? — >■ leads 
to kinetics of pointlike particles or fragmentation kinetics [25l [26l [30]. Though we 
limit ourselves to size ratios R > 1, the pointlike kinetics corresponds to the limit of 
R ^ oo, followed by a rescaling of the linear dimensions by the length of the large 
particle. 

To properly characterize the jammed state, we start by introducing some 
definitions and relations. The basic quantity to characterize is the probability 
distribution of the interparticle distance defined as 

„ , , , Number of empty intervals of size x within ]x,x + dx\ , , , 

Pti(x)dx = ^ / : -dx. (1) 

iotal number of empty intervals 

Thus, Pij{x)dx at the jammed state has the property 
1 

P^{x)dx = 1. (2) 

Discriminating all gap types in terms of pairs of consecutively adsorbed segments 
provides four distinct cases, namely AA, AB, BA, and BB. Density distribution 
functions Pb{x), with b € {AA, AB, BA, BB}, i.e., for each of the gap types, are 
defined, from which the following relation is obeyed 

P0(x) = Paa{x) + 2Pab{x) + Pbb{x), (3) 

where the equality Pab{x) = Pba{x) was used. Note from equations ([2]) and ([3]) 
that the Pb{x) are not normalized. Keeping in mind the above definitions, it is now 
straightforward to reckon moments of any order of the gap-size distribution functions, 
defined by 

= ^"f"^-^"^^" , (4) 
/o-Pa(a;)dx 

with a G {(l>, AA, AB, BA, BB} and n = 0, 1, 2, . . . , though in the present work, 
we restrict ourselves to moments with n < A. Using equation ([3]) one can relate the 
moments of the gap-size distributions with the corresponding moment of the global 
gap-size distribution function given by equation ^ yielding 

= (x")^^ f PAA{x)dx 



)ab / PAB{x)dx 

Jo 



^BB / PBB{x)dx, (5) 
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again, defined for all values of n = 0, 1, 2, We also compute the cumulants, kJ^, 

of the gap-size distributions defined as, 

lnG„(fc) = f;i^C, (6) 



m 

m—l 



where Ga{k) is the so cahed characteristic function defined by [JT], 

Jo Pa[X)dx 

Therefore, from equations (|4]), and ([7]), one derives the first four cumulants as 

< - (^)., (8) 
^^, = {x\-{x)l (9) 
4 = (x3)„-3(a;2)„(x)„ + 2(x)3, (10) 

< = (x^), - 4(:r3),(x), - 3{xY^ + 12(x2),(x)2 - 6(x)t (11) 

where Ki is just the mean value and K2 is the variance. While the first and 
second cumulants have well-established meanings, those of the third and fourth order 
cumulants are put in terms of the skewness. 



and kurtosis. 



Ka = -^, (13) 

respectively |31| . in order to facilitate their interpretation. For later convenience, we 
also take the opportunity to define the dispersion, aa = ^/k^, used in the analysis of 
the results provided in Sec. [3] 

2. The gap-size distribution functions 

Since the probability that a segment is deposited in a given gap is proportional to 
the effective size of the gap, the probability distribution describing the interparticle 
distance should be reasonably approximated by an exponential distribution. The 
effective gap size is defined as the difference between the actual size of the gap and 
the size of the segment attempting adsorption, so it represents the actual length 
available for deposition. In the jamming state, the state where no more particles 
can be deposited, the probability of having a gap of length larger than one is zero. So, 
we propose as the ansatz for the probability distribution function of the interparticle 
distance in the jammed state, an exponential function truncated to the interval [0, 1[, 

PJx,.{^"<-"-"'- lfl<\ (14) 

where aa is the single parameter used to characterize the distribution function of the 
interparticle distance. Since Pid{x) is normalized (eq. ([2|)), one has 

A^\a^) = = 1. (15) 

an, 
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Figure 1. Plot of Oa, with a G {AA, AB, BB}, as a function of size ratio 
computed from first moment of Monte Carlo simulations 1281 . Each of the 
AA-, AB-, and SB-gap types are represented by circles, squares, and triangles, 
respectively. 



Some considerations are due before proceeding to the derivation of the various 
cumulants. The idea of using the proposed ansatz (eq. (jl4p ) is to describe the jammed 
state. Though a truncated-exponential ansatz best describes uncorrelated deposition 
events, it allows us to assess the presence of correlated deposition events at the jammed 
state. The description of the full kinetics towards the jammed state would not be 
possible with the proposed ansatz. In particular, during the adsorption process, the 
upper limit of the integral in eq. ([2]) would have to be extended to infinity. Moreover, at 
early times, the kinetics of deposition is uncorrelated and an exponential distribution 
represents a better fit [6]. Thus, we are interested in the jammed state when systematic 
elimination, by the adsorption of particles, of all gaps of length a; > 1 have taken place. 

From eq. Q the various moments of the gap-size distributions in eq. (|14|) are 
calculated. The first moment equals the first-order cumulant (eq. ([8])), 

< = (^)a = -^ + — . (16) 



1 



The second moment is 



and the dispersion is 



aa{2 + Ua) 



1 



1 



(17) 



(18) 



y 2 — 2cosh(aa) 

To better characterize the distribution function, as previously referred to, one 
resorts to the calculation of higher moments. The third- and fourth-order cumulants 
are reckoned in terms of the moments and used for the definition of the skewness 
and kurtosis, which are quantities with a more intuitive meaning. The skewness 
yields information about the asymmetry of the distribution, while the kurtosis 
provides information about the peakness or flatness of the distribution function. For 
completeness, we provide the relations for the third and fourth moments, 

Q;a(6 + aa(3 + aa))' 



an 



6 + 



1 



(19) 
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Figure 2. Plots involving cumulants up to the fourth order of the gap-size 
distribution functions for each gap type as a function of the size ratio, according to 
equations l|16|l . I|18| l. I|23| l. and l|24p . a) Average and b) dispersion of the distance 
between pairs of segments, c) Skcwness. d) Kurtosis. The legend convention 
follows that presented in Fig. [l] 



and 



24 



aa(24 + a„(12 + aa(4 + aa)))' 



1 - e° 



(20) 



After some tedious, though straightforward, algebra we obtain the third and fourth 
cumulants, 

2 1 



and 



" a3 4 



6 1 

^ ~ 8 



coth I — I csch' 



(2 + coshK))csch4(^). 



(21) 



(22) 



From eqs. ^ and (HH), and using eqs. (HH]), and ([22]), we obtain the 

skewness 

4- i cosh {aa) + al coth {^f-) 



+ 2=2^EM^)i'^ + "a - 2 cosh («„)) ' 



(23) 



and the kurtosis 



J24 + al) cosh (gj + 2(^9 + - 3 cosh (2c^,)) 
(2 + a2„2cosh(a,))2 ' 

for the truncated exponential distribution function defined in eq. (jl4p . 
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Figure 3. Plots involving the ratio of cumulants of tiie gap-size distribution 
functions of the same order. The plots provide a direct comparison of the results 
from the truncated exponential ansatz to those from Monte Carlo simulations 
as explained in the text, a) Dispersion, b) Skewness. c) Kurtosis. The legend 
convention follows that of Fig. [T] 



3. Results and Discussion 

In a previous work |28| , we have performed an extensive Monte Carlo simulation study 
of the present random sequential model. Here, we utilize results from that work, 
specifically from the first moment as defined in eq. (|16p , to determine the values of a 
for all gap types and size ratios ranging from 1 to 20. Figure [1] shows the computed 
values of a as a function of the size ratio. Heretoforth, the values of the size ratio 
match those of the Monte Carlo simulations to simplify comparison [28] . For the BB 
gap type a monotonic decrease of a with the size ratio is observed. The AA and AB 
gap types are characterized by maxima for values of the size ratio at 1.30 ± 0.05 and 
1.20 ±0.05, respectively. 

To characterize the system, in Fig. [2] the average distance, dispersion, skewness, 
and kurtosis are computed, respectively, with eqs. (IT^ . (f^. and ([M)) as a 

function of the size ratio. For the AA and AB gap types, a non-monotonic behavior for 
all these quantities is observed. For the BB gap, the average distance and dispersion 
increase with the size ratio as shown in Fig. [DJa) and (b), while the skewness and 
kurtosis decrease monotonically as shown in parts (c) and (d) of the same figure. For 
the average distance and dispersion (Fig. [D^a) and (b)) there is a minimum for the 
AA gap at a value of the size ratio of 1.30 ± 0.05 and of 1.20 ± 0.05 for the AB gap. 
For the skewness and kurtosis (Fig. [2Jc) and (d)) there is a maximum for the AA gap 
at 1.30 ± 0.05 and at 1.20 ± 0.05 for AB gap. The values of the maxima and minima 
for the AA and AB gap types coincide with the minimum values of a as a function of 
the size ratio. 

We compare results using the truncated exponential ansatz with simulational ones 
found in Ref. [28]. All computed quantities possess the same qualitative behavior as 
that obtained from an extensive Monte Carlo study. The exception regards the values 
of the dispersion at sizes ratios near unit. In this situation, streaks of alternating A 
and B segments are common, since two consecutive BB segments might still have a 
gap of length larger than unit, where a single A segment can fit in. In fact, at size 
ratios lower than 1.55, the population of BB gaps remains significant, as compared 
with the A A gap one. However such streaks leading to gaps of size near unit, which we 
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denoted as snug fits in Ref. j28l , are not possible to be described by a simple truncated 
exponential as the latter does not account for influences due to the history of several 
segment deposition events. The ratios of the values of the cumulants of the same 
order from the truncated exponential ansatz to those from Monte Carlo simulations 
are show in Fig. [3l The BB gap shows better compliance with proposed ansatz. In 
a sense, this result is not entirely surprising, since the larger segments are deposited 
at early times. Consequently, there will be less interaction with previously deposited 
segments. At late times, i.e., in the asymptotic regime, the A segments will fill the 
gaps left out by the B segments, particularly as the size ratio increases. Improved 
compliance to the truncated exponential is also observed for the AA and AB gap types 
with increasing values of the size ratio for similar reasons to those of the BB gaps. As 
the size ratio increases, the B segments have a sharper decay in the rate of adsorption 
due to single A segments blocking their adsorption. Consecutive B segments also leave 
larger spacing between them, on average. Consequently, the adsorption of A segments 
follows more closely the argument in favor of a truncated exponential presented above. 
At size ratios close to one, the results provided by the ansatz deviate more, since the 
populations of the various gap types are nearly equilibrated and one expects the 
equations describing these gap types to be more coupled, i.e., to be more dependent 
on multisegment deposition events. 

4. Final remarks 

In short, we introduced a novel approach with potential for a semi-quantitative analysis 
of models where the interparticle distribution function of the distances is of relevance. 
Specifically, the method simply uses a single moment to obtain information for the 
other three. However, it requires some prior knowledge of what can constitute a 
good ansatz. We present a coherent study of the dependence on the size ratio up 
to the fourth moment. We propose a truncated exponential function for the gap- 
size distribution and compute the mean gap size, dispersion, skewness, and kurtosis. 
Qualitatively, the results for the proposed ansatz are in good agreement with the 
results obtained from simulational work. However, the dispersion of AB gap type at 
small values of size ratio, the exponential function does not reproduce the behavior of 
the distribution function, which can be explained by snug fit events. Finally, we point 
out that the present approach can be applied to a broader range of models involving 
the irreversible deposition of particles on a line. 
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